clc;
clear;
A=[-71,-65,-81,-46;
    75,89,117,50;
    0,4,8,4;
    -67,-121,-173,-58];
D=eig(A),
[v,d]=eig(sym(A))